New Precise Measurement of the Pion Weak Form Factors in 7r+ — > e+i/7 Decay 
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We have measured the tt"*" e^w) branching ratio over a wide region of phase space, based on 
a total of 65,460 events acquired using the PIBETA detector. Minimum-x^ fits to the measured 
{E^+,E^) energy distributions result in the weak form factor value of Fa = 0.0119(1) with a fixed 
value of Fv = 0.0259. An unconstrained fit yields Fy = 0.0258(17) and Fa = 0.0117(17). In 
addition, we have measured a = 0.10(6) for the dependence of Fv on g^, the e'^iy pair invariant 
mass squared, parametrized as Fv{(f) = -Fv(0)(l + a ■ q^). The branching ratio for the kinematic 
region E^ > lOMeV and 6^+^ > 40° is measured to be B^^^ = 73.86(54) x 10"*. Earlier deviations 
we reported in the high- i?^. /low- kinematic region are resolved, and we find full compatibility 
with CVC and standard V—A calculations without a tensor term. We also derive new values for the 
pion polarizability, as = 2.78(10) x 10^* fm'^, and neutral pion lifetime, r^o = (8.5 ± 1.1) x 10^^^ s. 



PACS numbers: 14.40. Aq, 13.20. Cz, ll.30.Rd 

Radiative pion decay 7r+ e+j/7 (also denoted ^02^), 
where 7 is a real photon, offers the best means to study 
pion form factors at zero momentum transfer. Pion form 
factors provide critical input to low-energy effective theo- 
ries of the strong interaction, such as chiral perturbation 
theory (ChPT). Previous experimental studies of the 7ro27 
decay were constrained by relatively low event statistics 
and limited kinematic coverage, leaving room for specu- 
lative interpretations including an anomalously large ten- 
sor interaction term. In this Letter we report precise new 
results of a combined analysis of data recorded with the 
PIBETA detector P, removing the limitations and am- 
biguities found in previous work. 

In the standard description, the 7ro2'y decay ampli- 
tude consists of the inner bremsstrahlung (IB), structure- 
dependent (SD*), and interference terms. The SD terms 
are parametrized by vector {Fy) and axial vector {Fa) 
form factors (see Ref. [2] for a review). The conserved 
vector current (CVC) hypothesis relates Fy to the 7r° 
lifetime [lll[5], yielding Fy = 0.0259(9) [6J. ChPT cal- 
culations [3 [HI [HI OLQJ give a value for Fa in the range 
0.010-0.012. 

Due to the above noted limitations, authors of previ- 
ous experimental studies fixed Fy at the CVC-predicted 
value and evaluated Fa, with resulting relative uncertain- 
ties ranging from ±12% to ±56% [HJ [H [13 El [13 [H] . 
None of these measurements was sensitive to the form fac- 
tor dependence on the e'^v pair invariant mass. Because 



of inconsistencies in the available data, a relatively large 
non-{V—A) tensor contribution to the SD terms was con- 
sidered (see Ref. [17] and references therein). Our first 
measurement |18j reported a significant improvement in 
the accuracy of Fa, but also noted a substantial deficit 
of observed events in the high-i?..y/low-£^e+ kinematic re- 
gion, still leaving room for non-iV—A) admixtures. 

We performed measurements at the irEl beam line 
at the Paul Scherrer Institute (PSI), Switzerland, using 
a stopped tt"*" beam and the PIBETA large-acceptance 
Csl calorimeter with central tracking [T. A total of 
over 2.0 x 10^'^ tt+'s stopped in the active target dur- 
ing the two measuring periods were used in the analysis. 
The first data set (1999-2001) was collected during the 
TT"*" —> TT^e'^v (ttis) branching ratio measurement (19] with 

8 x 10^ pion stops per second. A second sample was 
recorded in 2004 with an eight times lower stopping rate. 
Details of the detector design, performance, and our ex- 
perimental methods are presented in Refs. [T1 IT51 [TOll^n] . 

The PIBETA detector measured shower energies, 
and directions of the positron and photon, thus over- 
determining the kinematics of the final three-body state. 
Key for the 7102^ measurement were the one-arm and two- 
arm high-threshold triggers, IHT and 2HT, respectively, 
requiring that at least one (IHT) or two (2HT) showers 
register above an effective high threshold of ~ 54MeV 
deposited in the calorimeter, thus reducing a copious e'^ 
background from the ir-fi-e decay chain with a continuous 
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TABLE I: Total energy ranges for kinematic regions used in 
measurement {I-III, in measured MeVexp) and branching 
ratio evaluation {A~-0, in physical MeV). For regions B, C, 
and O, the e"''-7 relative angle is restricted to 9^+^ > 40°; this 
condition is automatically satisfied for regions I-III and A. 
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energy spectrum ending at ~ 52.8 MeV. In the 99-01 data 
set the IHT trigger was prescaled 16- or 8-fold. Our 7102^ 
data set covers three non-overlapping kinematic regions 
given in Tab. [ij region / data were acquired with the 
2HT trigger, while the IHT trigger mapped out regions 
// and ///. In all three regions the 71027 decay is strongly 
dominated by the SD+ term proportional to {Fy -I- F^)^ 
[5]. Candidate TTc2-y events exhibit one neutral shower in 
the calorimeter coincident in time with a positron track. 
For the events with more than one available combination, 
i.e., more than one neutral shower or positron track, the 
pair most nearly coincident in time was chosen. Events 
incompatible with the 7rc27 kinematics were rejected in 
the final data sample. 

Due to the nature of the detector resolution function, 
the experimental kinematic region contains a fraction of 
events shifted from a broader kinematic range. We there- 
fore used the reconstructed events in regions /, //, and 
/// to evaluate the branching ratios B'!^^^ for the larger 
regions A, B, and C, respectively, and O, the combina- 
tion of all three regions, by means of Monte Carlo (MC) 
extrapolation. On the other hand, we calculated the the- 
oretical branching fractions B^l?2^ for the same regions 
numerically using only the input parameters {Fy, Fa, a), 
as described below, and the standard description of the 
decay, including radiative corrections [21]. 

Accidental background was dominated by positrons 
from 7rc2 and 7r-/i-e decays accompanied by an unrelated 
neutral shower. It was corrected for by subtracting his- 
tograms of observables projected using the out-of-time 
cut Aiout '■ (5ns < \t^+ — t^l < 10ns) from the in-time 
histograms projected via the cut Atiii : {\tf,+ —ty \ < 5 ns). 
The 2004 data sample with lower stopped beam intensity 
provided an improvement of the signal to accidental back- 
ground ratios from 1.7:1 to 10:1 in region // and 4.4:1 to 
30:1 in region ///. 

Non-accidental background sources were predomi- 
nately TTfj events for which one tt" decay photon converts 
in the target, producing a charged track in the detector. 
Having measured the net yield of the decay events in 
our detector, we have used a MC simulation to calculate 
the fraction of tt^ decay events misidentified as 7rc2-Y. This 
class of background events contributes ~ 6.0% of the 



signal in the kinematic region / and negligible amounts 
in regions // and ///. The tt^ contamination was sub- 
tracted in the calculation of 77027 yields. 

In order to evaluate the branching ratio of the 77027 
decay in any given kinematic region, for normalization 
we used the 7ro2 events, recorded in parallel via the one- 
arm HT calorimeter trigger. This procedure assures the 
cancellation of the imprecisely known total number of 
stopped 7r+'s, and Ctrac; the combined tracking efficiency 
of e"'"'s. No statistically significant e"*" energy depen- 
dence was observed for ef^^^ I 



branching ratio for the tt^ 
lated from the expression 



2]. Thus, the experimental 
e'^v"/ decay can be calcu- 
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where A'^027 (-^7rc2) is the number of the detected tt, 
{tTc2) events, ^7ro27 (^7rc2) is the experimental acceptance 
for the given decay type incorporating the appropriate 
cuts, and i3,re2 = 1.2352(1) x 10"'' is the theoretical 7ro2 
branching ratio p3l l24j . The factor p^o2 denotes the 
prescaling of the one-arm HT trigger; it applied to the 
1999-2001 data set only. 

The experimental acceptances were calculated in a 
GEANT3 based MC program 25J. This simulation incor- 
porates a precise description of the PIBETA detector ge- 
ometry as well as the electronics response of the experi- 
mental setup. Figure [1] demonstrates the match between 
the MC simulation and the background-subtracted data 
in regions /, // and ///. The number A'^02 was extracted 
from an independent fit of the time distribution of the 7ro2 
event candidates. Upon imposing appropriate kinematic 
cuts and performing the background subtraction proce- 
dure we have reconstructed 35,948 ± 194 (0.54%) events 
in region /, 16,246 ±331 (2.0%) events in region // and 
13,263 ± 161 (1.2%) events in region ///, where num- 
bers in parentheses are fractional statistical uncertainties 
obtained after accidental background subtraction. The 
2004 data account for only 8 % of the total for region /, 
and about 43 % in regions // and ///. However, because 
of higher accidental background subtraction, the 2004 
events in // and /// are statistically more significant 
than the 1999-2001 events. More 1999-2001 events are 
included here than in the analysis reported in Ref. ^8] 
thanks to a broadening of the acceptance cuts made pos- 
sible by refinements in the treatment of the positron en- 
ergy deposited in the target. The acceptance cuts were 
adjusted so as to achieve minimal overall statistical un- 
certainty after background subtraction. The dominant 
sources of the systematic uncertainties are region depen- 
dent and are summarized in Table UT] Total statistical 
and systematic uncertainties for each region are summed 



in quadrature; they are given in Table III 

In order to extract the values of Fy, Fa, and a, and 
calculate the values of the branching ratios, we mini- 
mized the sum of differences between the experimen- 
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FIG. 1: Data points: background-subtracted 7r+ — > e dis- 
tribution of the kinematic variable A = {2Ec/m.^) s\v?{Oe^ /2) 
for regions / and II (top panel), and of the photon energy 
for regions / and /// (bottom panel). Solid lines: results of 
GEANT3 calculations for the best-fit values of Fv, Fa, and a. 



tal and calculated differential branching ratios over all 
phase space, taking into account the form factor depen- 
dence on the e^v pair invariant mass such that Fv{q^) = 
Fvmi+a-q^), FA{q^) = ^^(0), and = l-{2E^/m^), 
following the prescription of Refs. [H 13 HH]) where a is 
the slope parameter. The experimental branching ratios 
acquire a form factor dependence through the ac- 
ceptances At^c2i which are recalculated at every iteration 
step of the minimization. 

Both theoretical and experimental values of the 



TABLE II: Dominant sources of systematic and normalization 
uncertainties. They are: precise determination and simula- 
tion of the trigger thresholds, energy calibration for different 
types of particles, radiative corrections to the acceptance and 
number of the 7ro2 events. Remaining sources include uncer- 
tainties in the signal separation from prompt hadronic and nfi 
backgrounds. 
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0.20 
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0.13 


0.13 


0.13 


All other sources combined 


0.19 


0.02 
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Total 


0.79 


0.58 


1.38 



branching ratios are proportional to i37rc2, therefore our 
sum assures that the values of the extracted parame- 
ters are independent of the value of B.,,c2- 

For comparison with previous work, and to illustrate 
the sensitivity of our data, we performed a fit with fixed 
values of Fy = 0.0259 H, a = 0.041 [26^, yielding 
Fa ~ (119 ± 1) X 10~^ for the free parameter, with 
X^/n.d.f. — 0.96. This result represents a sixteenfold 
improvement in precision over the pre-PIBETA world 
average [27j. Table III and Fig. [T| show the excellent 



agreement between the measured integral and differen- 
tial branching ratios, and those calculated for the best-fit 
parameter values, refiecting the exceptionally good de- 
scription of our data by CVC and standard V—A theory 
in all regions of phase space. The previously reported 
anomalous shortage of events in region B [18 is thus re- 
solved. In this work the calorimeter energy calibration 
(CEC) for 7-induced showers was uncoupled from that of 
the e+-induced showers, resulting in a ^2 % difference. 
The clean 2004 data set revealed this discrepancy for the 
first time in the IHT tt^ decay event sample. Due to 
the limited statistical significance of such events in the 
1999-2001 set, it was not feasible to evaluate a separate 
7 CEC in our original analysis. Because of the extreme 
sensitivity of region // to the 7 energy scale, the effect 
was significant only in this region. More details can be 
found in Ref. |28j and in a forthcoming publication. 

Because of the strong dominance by the SD"^ term, 
our data are most sensitive to the combination Fa + Fy ■ 
This sensitivity is illustrated in Fig. |2] produced by a 
two-parameter {FA,Fy) fit. The best-fit value of Fa 
has an empirical linear dependence on the value of Fy 
parametrized as 

Fa = (-1.0286 • Fy + 0.03853) ± 0.00014 . (2) 

This expression will remain applicable as the CVC- 
predicted value of Fy changes thanks to future improved 
measurements of r^O) the neutral pion lifetime |29] . 

We also examined the question of uoyi-{V — A) contri- 
butions to the 7re27 decay by introducing tensor couphng 
following the prescription of Ref. [T7]. The optimal fit 
yields a single tensor form factor. Ft = (—0.6 ±2.8) x 



TABLE III: Best-fit 11^2-^ branching ratios obtained with 
Fv = 0.0259 (fixed), a = 0.041 (fixed) and Fa = 0.0119(1) 
(fit). Experimental uncertainties in parentheses refiect both 
the statistical and systematic uncertainties, given above. The- 
oretical uncertainties are dominated by the fit uncertainty in 
Fa. 
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B 
C 
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2.614(21) 
14.46(22) 
37.69(46) 
73.86(54) 



-Bthc X 10^ 
2.599(11) 
14.45(2) 
37.49(3) 
74.11(3) 



1.006(9) 
1.001(15) 
1.005(12) 
0.997(7) 
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FIG. 2: Contour plot of loci of constant for the minimum 
value xo plus 1, 2, and 4 units, respectively, in the Fa-Fv 
parameter plane, keeping the parameter a = 0.041. The range 
of the CVC prediction, Fv = 0.0259(9) is indicated. 



10-* , or -5.2 X lO-"* < Ft < 4.0 x lO^-* with 90% 
confidence. Ovir limits on Ft are in excellent agree- 
ment with Refs. j26l 130] , but are two orders of magni- 
tude more constraining than the results reported by the 
ISTRA group [ISIET]. 

Finally, in an unconstrained V—A fit we released all 
three parameters Fa, Fy, and a simultaneously, and 
obtained Fa = 0.0117(17), Fy = 0.0258(17) and a = 
0.10(6). These results are nearly identical to those of 
the two-parameter {Fa,Fv) fit presented in Fig. [2] and 
Eq. (|2|. 

In summary, our results: (i) are in excellent agreement 
with the CVC hypothesis prediction for Fy, (ii) repre- 
sent a fourteenfold improvement in the precision of Fy 
and a similar improvement in Fa, and (iii) provide the 
first ever measurement of the charged pion form factor 
slope parameter a, also found to be consistent with CVC. 
The PEN experiment [32J, currently under way at PSl, 
will add substantially to our TTc2-y data set. Our best-fit 
value of Fa agrees well with ChPT calculations, tending 
to the top of the reported range [IlIHlll]. However, a more 
precise measurement of Ttto is needed in order that the 
sensitivity of our data, expressed in Eq. ([2|, be put to full 
use in determining Fa- We use our form factor results 
to evaluate the pion polarizability and the ChPT param- 
eter sum Lg + as follows: using the one-parameter 
fit we obtain ue = -Pm = 2.78(2)cxp(10)fv x lO^-^fm^, 
and Lg -t- L^^q = 0.00145(1)oxp(5)fvi where the first un- 
certainty comes from the fit and the second from the 



current CVC-derived value of Fy- Alternatively, we get 
UE = 2.7 (^^) X lO-^fm^ and LI + L\a = 0.0014(1^) 
based on our unconstrained fit of Fa and Fy. In addi- 
tion, we use the latter fit result and CVC to make an 
independent determination of the neutral pion lifetime: 
tm = (8.5 ± 1.1) X 10-"s. 
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